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We consider magnetic impurities in a two dimensional superfiuid Fermi gas in the presence ol 
spin-orbit coupling. By using the methods of t-matrix and Green's function, we find spin-orbit 
coupling has some dramatic impacts on the effects of magnetic impurities. For the single impurity 
problem, the number of bound states localized around the magnetic impurity is doubled. For the 
finite concentration n of impurities, the energy gap is reduced and the density of states in the gapless 
region is greatly modified. 
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I. INTRODUCTION 

In the last few years, topological insulator (TI) and 
topological superconductor (TSC) [1, 2] have arisen a 
lasting interest in condensed matter systems. In such 
systems, spin-orbit coupling plays a fundamental role to 
make the TI and TSC different from the traditional in- 
sulator and traditional BCS superconductor significantly. 
In cold atomic field, the superfiuid of Fermi gases have 
been realized in 40 K [3] and 6 Li [4-6], furthermore, spin- 
orbit coupling has also been realized by the more recent 
experimental achievement of synthetic gauge field [ , ] . 
By Introducing the effect of spin-orbit coupling, new 
phases analog to TSC and new phenomena may emerge 
in cold superfiuid Fermi gases, e.g. topological superfiuid 
and topological phase-separation [9]. Since a cold atomic 
system is more adjustable, it can provide a more control- 
lable way to study such phases under different conditions, 
e.g. magnetic impurities. 

From TI and TSC, we know that spin-orbit coupling 
combining with Zecman field greatly change the energy 
spectrum and make topological phases emerge. There- 
fore, to figure out how the spin-orbit coupling affects the 
effect of magnetic impurity, which likes a local magnetic 
moment, is a very interesting problem. Recently, peo- 
ple have paid attention to studying on the Kondo effect 
in the presence of spin-orbit coupling [10-12], however, 
in this paper, we will study the effects of classical mag- 
netic impurity in superfiuid Fermi gases in the presence of 
spin-orbit coupling. Since cold Fermi gas system is much 
cleaner than metallic superconductor and the strength of 
pairing interaction can be tuned by adjusting the thresh- 
old energy of Feshbach resonance, obviously, it's a better 
platform to detect the effects of dilute magnetic impuri- 
ties. 

The effects of magnetic impurities are particularly im- 
portant in the case of superconductor, and so superflu- 



idity. Magnetic impurities which can flip the spins of 
electrons interfere with the pairing and reduce the gap. 
A sufficient concentration of magnetic impurities leads 
to gapless superconductivity, where the superconductive 
order parameter still exist but the BCS excitation gap is 
absent [13, 14]. Such a pairing-breaking effect is absent 
in the case of nonmagnetic impurities (known as Ander- 
son's theorem), so in the following, we only consider the 
effects of magnetic impurities. 

In this paper, we discuss the effects of magnetic im- 
purities in a two dimensional superfiuid Fermi gas in the 
presence of spin-orbit coupling. So far, nonmagnetic im- 
purities have been realized in Bose gases by using other 
species of atoms [15] and laser light [16]. Recently, a 
single spin-J, placed in a Fermi sea of spin-f atoms is 
observed to behave as a mobile impurity to form Fermi 
polaron [17, 18]. To realize magnetic impurities in a su- 
perfiuid Fermi gas, Y. Ohashi [19] suggests and confirms 
that nonmagnetic impurities accompanied by localized 
spin-f can be viewed as magnetic impurities. Based on 
these, we believe it is not far away to detect the effects of 
magnetic impurities in a superfiuid Fermi gas with spin- 
orbit coupling in real experiments. 

The paper is organized as follows. In Sec. II, we con- 
sider the case of single magnetic impurity and calculate 
the energy of bound state in the presence of spin-orbit 
coupling. In Sec. Ill, we generalize the single impurity 
case to the one of concentration n. In Sec. IV, a conclu- 
sion is given. 



II. SINGLE IMPURITY MODEL 
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Now, we consider a two dimensional ultra-cold atomic 
system with single magnetic impurity and spin-orbit cou- 



pling, which is described by 
H = Hq + Hj, 



'imp: 



= J2 £kC kv Ck ° + 

key 



E 

k 



a R k ( 



where p(Q) is the density of states at the Fermi level in 
the normal state. Based on Eqs.(5) and (7), the poles of 
t(w) arc given by (see details in Appendix) 
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where e k = k 2 /2m — /x, is the strength of spin-orbit 
coupling, c ka (c kcr ) is creation(annihilation) operator with 
pseudo-spin er, A is the order parameter, U is the on-site 
interaction strength and s z = ±s is the number of excess 
spin- up atoms localized around impurity (see Refs.[19]). 
For simplicity, here we completely neglect the effect of 
the local spatial variation of the order parameter around 
a paramagnetic impurity. 

In the Nambu spinor presentation ^1 = 



(4t' c H' c -H> c - 



Hn can be rewritten as 



Ho = -5>iff(fc)*fc, 

fe 

H(k) = e fe T 3 - ApiT 3 + a R k 



,(2) 



where pi and Ti(i = 1,2, 3) are Pauli matrices in particle- 
hole space and spin space, respectively, and r± = (n ± 
ir 2 )/2. According to Ref.[20] or the method of the co- 
herent state path integral [21], we can have 

Go(k.ui) = t —r-. : r, 

u - e k r 3 + ApiT 3 - a R k{e l ^ k T + + e^^ fc r_) 



where A = ([/s7rp(0)/2) 2 . From Eq.(8), it's easy to see, 
no matter how small A may be, two bound states always 
appear, this result is quite different from the case in the 
absence of spin-orbit coupling. However, for a R = 0, 
there are only two poles in the energy gap, and the poles 
reduce tow = ±[(1 — A)/{1 + A)]A, just the same result 
as in Rcfs.[19, 20, 22]. Increasing the strength of spin- 
orbit coupling, the energy of the low-lying state which is 
localized around the impurity goes to zero. From Eq.(8), 
we can see, if we increase the strength of U and keep 
A < 1, the combined effects of magnetic impurity and 
spin-orbit coupling is strengthened, and the energy of 
the bound state will go to zero faster than the case in 
the absence of spin-orbit coupling. 



III. IMPURITY OF CONCENTRATION N 

Now, we generalize single impurity to the case for fi- 
nite concentration n. After averaging over the random 
distribution of the impurities, the Green's function re- 
covers the translational invariance, and the renormalized 
Green's function is given by 



(3) G(k,u) 
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Go(fc,w) is related to the full Green function G(k,u>) in 
the following form 

G(k,k',u) = G Q {k,Lo)8 kk , +G {k,u)t{io)G (k\io). (4) 

After averaging over the direction of the spin, the non-flip 
part of the t- matrix, t{uj), is given by 
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where Q and A are the renormalized frequency and or- 
der parameter, respectively. u> and A are determined by 
Dyson equation 

G(k,u) = G (k,u) + G (k,uj)Z(k,u)G(k,u), (10) 

another useful form of this formula is 

Gi^uj)- 1 = G (k,uj)- 1 -Z(k,oj), (11) 

where Go(k,u), the Green's function in the absence of 
impurities, is given by Eq.(3). In Born approximation, 
the self energy Y,(k,u) is given by [14] 



We are interested in the localized excited state in the 
energy gap A, so we must find poles of t(ui). For \u\ < A, 
by integrating out k in Eq.(6) (see details in Appendix), 
we obtain 
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where t s is the relaxation time and is given by 
t s = Trp(0)ns 2 U 2 . 
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By substituting Eq.(12) into Eq.(ll), we obtain 
U)i = uji + 
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where Ai = A 2 = A. They arc quite different from the 
case without spin-orbit coupling, as before, with spin- 
orbit coupling, the number of poles of t-matrix has in- 
creased from two to four. By introducing a new auxiliary 
parameter u, defined by Ui = Wj/Aj (i = 1,2), Eqs.(14)- 
(15) are reduced as 
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Plotting the above equation in the uw-plane (as shown 
in Fig.l), we can see the case in the absence of spin-orbit 
coupling is quite different from the one with spin-orbit 
coupling. 

First, for £ < 2 case, with the spin-orbit coupling, 
there are two poles in Eq.(16), one is at u c = yj\ — C 2 /4, 
and the other is the original one at u c = 1. In the ab- 
sence of spin-orbit coupling, for v = 0, uj goes through a 
maximum and afterwards decrease monotonically, this is 
also the characteristic of Fig.l(a)(c)(d). However, there 
is no such characteristic in Fig. 1(b), after uj goes through 
a finite maximum, uj will not decreases monotonically, in 
the opposite, when u goes to u c , uj is positive and di- 
verges. So only the maximum value of uj in Fig. 1(a) has 
the meaning of the energy gap w g (in the following, we 
only consider this case). It is determined by the condition 
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This equation is difficult to solve. However, by numerical 
calculating, we find that if v be sufficiently small, like 
v = 0.05 in Fig.l, we can safely ignore the last term in 
Eq.(17) (see Fig.2), and obtain 
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FIG. 1: (Color online) The parameters of Eq.(16) are: £ = 
0.5. The values of Vi are shown in the figure, (a) There is a 
maximum for every fixed vi, and the value of the maximum 
decreases monotonically with vi. (b) There is no maximum, 
when it goes to the critical value u c , the value of uj/ A diverges. 

(c) When u £ (u c , 1), to /A decreases with u monotonically. 

(d) There is a maximum for every fixed non zero V2, and the 
value of the maximum decreases monotonically with V2 




U 



FIG. 2: (Color online) The parameters of Eq.(17) are: ( = 
0.5. f(u) = By increasing v from v = 0.05 to v = 0.15, 

wo, the it corresponds to the maximum of ^, changes slowly. 



We see that the solution of Eq.(17) exist only when £ < 1. 
In fact, even we consider the effect of spin-orbit coupling 
and modify u to (1 — a — C 2 / 3 ) 1 / 2 , to first order, the 
expression of uj g in Eq.(16) will not change. Based on 
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above, we can expand Eq.(16) around uq and obtain 
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Second, for 1 < ( < 2 case, lj decreases monotoni- 
cally as u increases between [0,u c ). This corresponds to 
Lu g = 0. Because of the presence of spin-orbit coupling, 
from Eq.(18), we can see that, when £ < 1, LJ g can reach 
zero. This indicates that spin-orbit coupling modifies 
strongly the critical concentration n c ( n c corresponds to 
£ = 1 in the absence of spin-orbit coupling), where the 
gap in the energy spectrum uj g vanishes and the gapless 
superfluidity appears. 

For small values of w, an asymptotic expression of u 
based on Eq.(15) is given by 



u = iV(C + a) 2 - 1 
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where a is a small constant determined by 
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For sufficiently small v, a can be ignored, and Eq.(20) 
can be reduced as 
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In the following, in order to show what are changed 
by the spin-orbit coupling, we calculate the density of 
states, which can be measured by STM. The density of 
states is given in terms of the Green's function by 
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Making use of Eq.(16) (here we can also safely ignore the 
last term), the above formula reduces to a more conve- 
nient form: 

N s (lo) = iV(0)C _1 Imu. (24) 
First, let us consider £ < 1 case, by using Eq.(19), we 
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FIG. 3: (Color online) Density of states plotted for different 
values of G The parameter v is set to 0.1. 



substitute the expression of u into Eq.(24) and obtain 
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From Eq.(25), we see that oj g gives the threshold fre- 
quency of the density of states, and spin-orbit coupling 
reduces the value of Ns(lj) for lj > uj g . This is obviously 
right, since from Eq.(18), we see spin-orbit coupling re- 
duces the value of ui g . To guarantee J N s (u))cLu invariant, 
the reduction of Ns (lj) is necessary. 

Second, for ( > 1 case, by Eq.(22) we obtain 
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here we have reconsidered the second order of lj/A. From 
Eq.(26), we can see N s is finite at lj = and given by 
N S (Q) = (1 - C 2 ) 1/2 N(0). This is just a characteris- 
tic of the gapless region, where the energy spectrum oj g 
vanishes though the gap parameter A is not zero. This 
is an example of gapless superfluidity, the existence of 
Cooper pairs without the existence of an energy gap in 
the spectrum of excitations. 

From Eq.(26), we also find that, for lj/A sufficiently 
small, Ns{lj) oc lj, which is quite different from Ns(lj) oc 
lj 2 in the absence of spin-orbit coupling [14]. From all 
of the above, it's obviously that the spin-orbit coupling 
affects the effects of magnetic impurities a lot. 
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IV. CONCLUSIONS 

In this paper, we have investigated the effects of mag- 
netic impurities in a two-dimensional supcrfluid Fermi 
gas with spin-orbit coupling. In the presence of spin-orbit 
coupling, by using the methods of t-matrix, we find that 
the number of bound states has doubled comparing to 
the case without spin-orbit coupling. By using the meth- 
ods of Green's function, we also find that the spin-orbit 
coupling changes the energy gap in the energy spectrum 
and density of states. In the gap region, the larger the 
concentration, the stronger the impact is given by the 
spin-orbit coupling. In the gapless region (here limited 
to £ > 1), we find spin-orbit coupling changes the rela- 
tion between N(ui) and ui from N(ui) oc uj 2 to N(uS) oc u> 
in the small u limit. 

Since we only consider the spin-flip effects of the mag- 
netic impurities and totally ignore the spin-exchange ef- 
fects, the impurities are classical impurities and as a re- 
sult, the Kondo effect is absent. In future, such a system 
including spin-exchange term may be used to explore the 
effects of spin-orbit coupling to the Kondo problem in 
cold atomic systems. In addition, the effects of magnetic 
impurities in imbalance Fermi gas with spin-orbit cou- 
pling are also worth exploring. 



Acknowledgments 

We thank Liang Chen for helpful discussions. This 
work is supported by NSFC Grant No.10675108. 



Go(k,ui) is a 4 x 4 matrix, after we write down its entities, 
by using the formula 



(A2) 



we can integrate out k directly, and obtain the entities 
of F(ui). However, because the presence of spin-orbit 
coupling, the energy spectrum has greatly modified, and 
this makes the integration difficult. To avoid this dif- 
ficulty, we assume the strength of spin-orbit coupling, 
an << Vf, where Vf = kp/m is the Fermi veloc- 
ity. Following the assumption k — kp ~ me/kp made 
by G. Rickayzen when he solved the problem of impuri- 
ties in metal[23], assuming the particle- hole symmetry of 
Fermi band and approximating the normal-state density 
of states by the value p(0) at the Fermi level, we obtain 



F(u) = -7T/?(0) .„ = n -. (A3) 



VA 2 - ui 2 



Using this formula and Eq.(5), we obtain the equation 
for the bound-state energies as 



= dct[l- AF(cj) 2 ], 
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Appendix A: Appendix 



where A = (U sir p(0) / 2) 2 . By directly solving a eigen- 
value problem and assuming ma 2 R / A << 1, we obtain 



To obtain Eq.(7), we first write down the Green's func- 
tion in the absence of a magnetic impurity, 
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